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A relativistic theory of gravity has recently been proposed by Bekenstein, where gravity is me- 
diated by a tensor, a vector and a scalar field, thus called TeVeS. The theory aims at modifying 
gravity in such a way as to reproduce Milgrom's Modified Newtonian Dynamics (MOND) in the 
weak field, non-relativistic limit, which provides a framework to solve the missing mass problem 
in galaxies without invoking dark matter. In this paper I apply a covariant approach to formulate 
the cosmological equations for this theory, for both the background and linear perturbations. I 
derive the necessary perturbed equations for scalar, vector and tensor modes without adhering to a 
particular gauge. Special gauges are considered in the appendix. 
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I. INTRODUCTION 

Bekenstein has recently proposed a relativistic theory of gravity where gravity is mediated by a tensor, a vector 
and a scalar field, thus called TeVeS |1|, aiming at explaining the missing mass problem. 

The missing mass problem is the longest standing problem of modern cosmology. It spans a wide range of scales, 
from galaxies to the cosmic microwave background. The problem is easy to state : the observed mass coming from 
all visible matter at the scales of interest, cannot account for the Newtonian (or General Relativistic) gravitational 
force observed acting on the same objects. The problem has a long history and manifests as discrepancies in the 
rotation curves of galaxies, motions of clusters of galaxies, gravitational lensing, and the absence of strong damping 
of linear perturbations on very large scales, to name a few. 

One could imagine that this missing mass is composed of baryons in objects other than stars, for example jupiter 
size planets or brown dwarves, collectively called MACHOS, or baryonic dark matter. These objects cannot be seen 
because they do not emmit light of their own. However microlensing studies did not detect the abundunce needed 
for these objects to make up for the missing mass 3]. Moreover the abundances of elements predicted by Big Bang 
Nucleosynthesis (BBN) give a matter density far below the needed mass density 4] . 

The popular approach to solving the missing mass problem, is to posit a matter component which does not interact 
with electromagnetic radiation and therefore cannot be detected by observing photons at various frequencies. Even 
though it cannot be seen directly, its presence is evident from the pull of gravity. Thus one attributes the extra 
gravitational force observed, to a "dark matter" component whose abundance is required to greatly exceed the visible 
matter abundance. Dark matter candidates have been traditionally split into "hot dark matter" and "cold dark 
matter", although an intermediate possibility, namely "warm dark matter" is sometimes beingconsidered. 

The earliest possibility considered for a dark matter candidate was a massive neutrino 0,13 j since neutrinos are 
particles which are known to exist as well as being very weakly interacting. However, massive neutrinos cannot be 
the dominant form of the dark matter. If the dark matter is composed of massive neutrinos then their mass must 
be at most 30 — 70eV for reasonable values of the Hubble constant, if they are not to overdose the universe 0. On 
the other hand the Tremaine-Gunn inequality gives a lower bound on the neutrino mass if neutrinos are to be 
bounded gravitationally within some radius. For example for dwarf spheroidal galaxies, their mass should be greater 
than ^ 300 — 400eT^ which is well above the cosmologically allowed mass range. Finally the recent Mainz and Troisk 
experiments from tritium beta decay, combined with neutrino oscillation experiments give an upper limit for the 
neutrino mass of around 2.2 — 2.5eV jS| . Massive neutrinos are therefore ruled out as a dark matter candidate capable 
of solving the missing mass problem. 

Cold dark matter, is composed of very massive slowly moving and weakly interacting particles. A plethora of 
such particles generically arises in particle physics models beyond the standard model quite naturally. The list of 
canditates is very long, ranging from light particle s lioll . supersymmetric particles to Kalutza-Klein modes 
and many more exotic objects. This subject (see for a nice recent review) has been studied in great depth and 
has been shown to agree with observations to a very good degree. Still, the actual nature of cold dark matter is left 
to speculation at the present time. 

The alternative approach is to point the finger to the laws of motion, or the gravitational field itself. This path 
was initially followed by Milgrom who proposed that for accelerations smaller than some acceleration scale oq, 
gravity departs from Newtonian gravity (which is still valid for large accelerations), in such a way as to explain the 
flat rotation curves of galaxies. It was thus dubbed Modified Newtonian Dynamycs (MOND). Alternatively it was 
shown that one can cast a modified dynamics theory into a modified gravity theory which provides essentially the 
same phenomenology 'l5l , although it might be possible that the two are not entirely equivalent '16] . Other modified 
gravity theories beyond the MOND paradigm have also been proposed (17i | but I do not consider them further in this 
paper. 

While MOND was successful as a simple phenomenological model in describing rotation curves, it had other serious 
problems. Taken at face value, it violates conservation of momentum, energy and angular momentum |l4j . This 
though stems from the fact that it is not a theory but rather an empirical law. Bekenstein and Milgrom have found 
a non-relativistic self-consistent realization of MOND, based on an aquadratic Lagrangian, thus called AQUAL [l5l |. 
However being non-relativistic, the theory could not make clear predictions about cosmological scales, for example 
the cosmic microwave background (CMB) or the formation of linear structure (the relativistic AQUAL proposed in 
the same paper suffered from acausal propagation of small perturbations). 

Nevertheless, even in the absence of a consistent relativistic MOND theory, several authors have tried to squeeze a 
cosmology out of MOND. Some authors based their calculations on general relativity to conclude that MOND is not 
compatible with cosmological observations ^ ^^ot so robust conclusion; the correct conclusion would have been 
that a baryonic, dark matter free universe, evolving under Einstein gravity cannot fit the cosmological observations. 
In other words they have shown the missing mass problem on cosmological scales. Nevertheless, those studies showed 
that relativistic MOND theories would have to give spectra very similar to dark matter cosmologies, if they are to 
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be considered as serious competitors. Given that relativistic MOND theories such as TeVeS have genericahy more 
parameters, this would tend to favour dark matter cosmologies when tests such as bayesian evidence are performed, 
unless their likelihood is strongly peaked about a very small region in parameter space, which is still an open question. 

Other authors 19] used GR to make predictions for early universe cosmology (for example the CMB) where MOND 
effects were argued to be negligible and then used heuristic arguments to make predictions about the growth of 
linear/non- linear structure. As 203 have shown, their CMB predictions are indeed compatible with the robust TeVeS 
calculations for some range of parameters, but can be quite different in other cases [2ll |. The matter power spectrum 
predictions exhibit a similar behaviour: the baryonic oscillations which were thought to be a MOND prediction, can 
be absent in TeVeS, just like in dark matter cosmology. 

Lue and Starkmann assumed that Birkhoffs theorem would be satisfied in a relativistic MOND theory. If that 
is true one can then build a cosmology out of MOND, in the same way one can derive the Friedmann equation in the 
matter era out of Newtonian gravity. What they found was that the growth rate of perturbations in the matter era 
in this MOND-like cosmology is slower than in GR cosmology with dark matter. Therefore if MOND is to be the 
limit of a relativistic theory which can successfully fit cosmological observations, it MUST violate Birkhoff's theorem. 
Indeed, TeVeS theory does violate Birkhoff's theorem; not only are spherically symmetric solutions not necessarily 
static, even the static solutions are not unique ,23,] . The bottomline is that to make robust cosmological predictions, 
consistent relativistic MOND theories are needed. 

Relativistic MOND realizations were constructed to overcome the acausal propagation of relativistic AQUAL 
typicall y b ased on two metrics which are conformally related via a scalar field. Phase coupling gravitation 
(PCG) |2J,|2^|2g is one such example. However the PCG parameters that could provide good MOND phenomenol- 
ogy were ruled out by solar system tests. Moreover, just like relativistic AQUAL, PCG cannot provide the observed 
gravitational lensing from visible matter alone. Part of the problem is the conformal relation between the two 
metrics [23, ll^- One can generalized this relation to a disformal one (2^ by including an additive tensor in the 
transformation, not related to the two metrics, for example build out of the gradient of a scalar field. However it 
was soon realized that any generalized scalar-tensor gravitation theory, even with a disformal relation between the 
two metrics in the theory, would produce less bending of light than GR and thus could not be used as a basis for 
relativistic MOND 

Sanders stratified theory [sj] manages to solve the lensing problem. Instead of using just a scalar field to disformally 
relate the two metrics in the theory, a vector field is used. The vector field in Sanders stratified theory is however non- 
dynamical, which contradicts the spirit of general covariance. This was solved by Bekenstein, who made the Sanders 
field dynamical by including an action for it (which is a special case of Jacobson-Mattingly [S^] Einstein-Ether theory 
action). The resulting theory was called Tensor Vector Scalar (TeVeS) gravitational theory, and was shown to 
provide MOND and Newtonian limits in the weak field non-relativistic limit, and was devoid of acausal propagation 
of perturbations. 

Ultimately, as every physical theory, TeVeS has to face as trop hysical and cosmological observations on every scale. 
Some strong lensing studies have already been performed jSJ]. The theory was also confronted with observations 
from our own galaxy . Alternatively one can test the theory with large scale structure observations or the Cosmic 
Microwave Background, both of which require the propagation of linear perturbations about a Friedman-Lemaitre- 
Robertson-Walker(FLRW) background. A first study in this direction has already been carried out The subject 
of the present work, deals with the formulation of linear perturbation theory in TeVeS, generalizing the equations 
given in po| to include curved spatial hypcrsurfaces, multiple gauges and all perturbation modes (scalar, vector and 
tensor modes). This opens up a whole new series of observational tests that can be performed on TeVeS, involving 
the largest scales in the universe. 

This paper is organized as follows. In the first section I give a short overview of the TeVeS theory, focusing on the 
action and the field equations. At the same time I introduce a somewhat simpler notation than the original TeVeS 
paper Q which can simplify both the actions and the field equations. Relativistic fluids are then introduced which 
are of particular importance to cosmology. The section concludes with the transformation of connections between the 
two frames associated with the theory. 

In the second section, I lay down a covariant formulation of the FLRW cosmology. The covariant equations are 
derived on the basis of the symmetries of the FLRW spacetime. A short discussion of the effective Friedmann equation 
follows, focussing on the time variation of the effective gravitational coupling strength, and the definition of the relative 
fluid densities. A specific choice of a coordinate system relevant to calculations for Cosmic Microwave Backround 
(CMB) anisotropics and large scale structure power spectrum, is given in the end. 

The third section takes on linear perturbations about the FLRW cosmology. The relevant tensors are perturbed 
covariantly, without adhering to a particular gauge, or perturbation mode. Thus the final perturbed equations contain 
scalar, vector and tensor perturbations. The use of the covariant approach is important due to the non-uniqueness of 
connections, which depend on which metric is used. Since the transformation of connections is derived in a covariant 
fashion, the final perturbed equations can be derived unambiguously. This section also gives relations of the perturbed 
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metrics in the two frames. 

In the fourth section, I take the perturbed field equations derived in the third section, and split them into irreducible 
parts by seperating out the scalar, vector and tensor perturbation modes. The resulting equations are still in a gauge 
non-fixed form, which makes it straightforward to check that they are indeed gauge invariant as expected. 

Finally the paper concludes in section six, where a summary is given of the results. 

The reader will also find useful the appendices. In particular appendix-ogives the perturbed equations for scalar 
modes in three different gauges while appendix-El gives a lot of intermediate steps in the derivation of the perturbed 
Einstein tensor. These steps were not included in the main part of the paper, to make it more readable, but are very 
useful to have when following the calculations. 

Throu gh o ut the paper I will use a signature +2 metric, and the curvature conventions of Misner,Thorne & 
Wheeler |43|. Greek indices are abstract tensor indices with no respect to any coordinate system. When writing 
tensor components in a particular coordinate system, Latin indices are used for the spatial part of the tensor with a 
"hat" on the index, while is used for the temporal part of the tensor. I will also use units such that the speed of 
light, Planck's constant divided by 2n and Boltzmann's constant are all equal to unity. 

II. FUNDAMENTALS OF TEVES 
A. Preliminaries 

TeVeS theory is a bimetric theory where gravity is mediated by a tensor field g^n with associated metric-compatible 
connection and well defined inverse g^^ such that g^'^Qpi, = S^^, a timelike (dual) vector field ^4^ such that 
g'^'^Afj^Ai, = —1, and a scalar field cj). Matter is required to obey the weak equivalence principle, which means that 
there is a metric g^j^ with associated metric-compatible connection V^, universal to all matter fields, such that test 
particles follow its geodesies. The tensor field g^^^ will be called the Einstein- Hilbert frame metric (see below) while 
gi_iu the matter frame metric. 

The relation between the four above tensor fields (when the field equations are satisfied) is 

g^.u = e-2"^g^^ - 2 sinh(20)A,, (1) 

with inverse 

gt.^ ^ _^ 2 sinh(2</)) A^A'' (2) 

where A^" = g'^'^A^. 

B. The action principle 

The theory is based on an action 5, which splits as 5' = Sg + Sa + + Sm, where Sg,SA,S^ and Sm are the actions 
for g^^, vector field A^, scalar field and matter respectively. 

The action for g^i/, Af^ and (f> is most easily written in the Einstcin-Hilbert frame, and is such that Sg is of 
Einstein-Hilbert form 

where g and R are the determinant and scalar curvature of g^i, respectively and G is the bare gravitational constant. 
Due to the complicated nature of the equations, the numerical value of G will not be the measured value of Newton's 
constant as measured on Earth. The precice relation between them depends on the spherically symmetric solution 
which apart from depending on the arbitrary function V (see below), is also not unique [lll23l|. 
The action for the vector field Af^ is given by 

Sa = Jd'x^i Kb F^,F>^^ - 2\{ApA>^ + 1)] , (4) 

where F^i, = 2V[^Aj,], F^^" = g'^'^g"^ Fap, A is a Lagrange multiplier ensuring the timelike constraint on and Kb 
is a dimensionless constant. 
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The action for the scalar field cf) is given by 

1 



Sdi 



16nG 



(5) 



where fi is a, non-dynamical dimensionless scalar field, and V{fj,) is an arbitrary function. The arbitrary function V is 
related to Bekenstein's function F as 



(6) 



where /i — SttGit^, (Jb being Bekenstein's auxiliary scalar field. Note that I have absorbed Bekenstein's constant ks 
into my definition of the function V. 

The matter is coupled only to the matter frame metric g^^ and thus its action is of the form 

(7) 



for some generic collection of matter fields x^- 



C. The field equations 

Variation of the action w.r.t the matter frame metric gives the matter stress-energy tensor as usual by 

1 



(8) 



where T^i^ is the standard matter stress-energy tensor. 

Variation of the action with respect to the three gravitational fields gives the field equations in the Einstein- Hilbert 
frame. 

The field equations for g^i, are given by [isf 



where G^i, is the Einstein tensor of the tensors y^,y and S^i, are given by 

+Kb 



1 



and 



S^, - T^. + 2(1 - e-^^)A^T^^^A,^ 



respectively, where V' = 

The field equations for the vector field A^^ are 

where the current is given by 

= (1 - e-4^)A^Tv 

The Lagrange multiplier is not arbitrary but can be calculated by contracting (|12|l with A^ and is given by 

A = KbA^V^F'^" + tiA^A''V^(f>V,(f> ~ SttG A'^j^ 
Inserting the above equation in H12|l one gets alternative field equations for the vector field as 



(9) 



(10) 
(11) 

(12) 
(13) 
(14) 

(15) 
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which do not expHcitely include the Lagrange multipher. 
The field equation for the scalar field <j) is 

V^r^ = SirGJ (16) 

where 

= ^ (gf^ - A" A") V^cj) (17) 

and where the scalar source J is given by 

J = e-^''' [g"" + 2e-^'''A''A''] T^^ (18) 

Apart from the field equations above, TeVeS theory has two constraint equations which are as follows. The first 
constraint is nothing but the timelike constraint on the vector field, 

r'A^A^ = -1. (19) 

which is found by varying the action with respect to the Lagrange multiplier A. The second constraint, fixes the 
non-dynamical scalar field /i in terms of the other fields of the theory. It is given by 

(^M-_^M^-)V^0V,0^_y' (20) 

which is found by varying the action with respect to /i. The second constraint equation must be inverted to find /i 
as a function of g^", and 0. Therefore the arbitrary function V and its derivatives are nothing but functions of 
kinetic terms for </>, contracted with g^^^ and A^^ . 

D. Fluids 

The stress-energy tensor of a fluid with density p, pressure P, velocity u^^ = g^^u^ and shear S^^ is 

T^^u = {p + P)Uf^u,y + Pg^^ + S^j,. (21) 

The velocity vector field is normalized with respect to the matter frame metric as g^'^u^Uu — u^^u^^ = — 1 while the 
shear obeys the two conditions = and u'^'Sfj,^ = 0. 

Using pi(l one obtains the contribution of the fiuid stress-energy tensor to the generalized Einstein equations as 

Si,u = ip + P) [uf^u^ -f 2(1 - e-*"^) + P [g^^ -I- 4 sinh(20)A^ A^] 

+ 2(1 - e-^^)A^^^^^A,) (22) 

Similarly equation (|13|l gives the contribution of the fluid stress-energy tensor to the vector field equations as 

ju = (1 - e-4*) [[p + P)A^'u^, + A^S^,] + 2 sinh(20)PA, (23) 

Finally the contribution of the fluid stress-energy tensor to the scalar field equation is obtained from p8(l as 



P- p + 2e-^'^{p + P){A^'Uf,)^ + 2e~^'^Ai'A''T,f,^ 



(24) 



The fluid evolution equations are obtained as usual from V fjT^,j = 0, where T'^^^ — g'^PPp^. Using (|21|) and 
expanding, gives 

u^^u^Wf^ip + P) + {p + P)u,V^uf' + {p + P)u^'WpU, + V,P + V^E^, = (25) 
where = ij^^'Ep^ . Contracting (|25() with u"^ gives the energy " conservation" equation as 

u^'VpP +ip + P)V^u'' = 0, (26) 
and subtracting (|26|) from (|25|l yields the momentum transfer equation as 

(p + P)u^y^u, + ((5^, + u.w^) V^P + V^E^, = 0. (27) 
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E. Transformation of connections 

I conlude this section by considering the transformation of connections from to and vice versa. This will 
come out to be very useful below, particularly in linear perturbation theory. 

Consider two metrics (7^1, and g^i, on a manifold M, with connections and respectively (not necessarily 
the two metrics of the TeVeS theory). The connections are required to agree on scalars, i.e. V^/ = V^/ for any 
/ € C°°{M). Acting on any form G T*M, the connections are related by 

V^u^ = V^M^ - D^^^ux (28) 

where the connection tensor -D^^^ is given by 

= ig^" (V^5p, + V,.9p^ - Vp5p,) (29) 
Using the metric relations Q and Q in H29|) one gets for TeVeS 

= 2J"(^V,)<^- [5p. + 2e2*ApA,]g"'5V/30 + 4A"A(^V,)0 

+ (1 - e^4^)A"V(^A,) + (e^^ - l)A^,,F^^" + 4smh\2^)A^^,F^f' ApA'^ (30) 

whereas with respect to V^, 

I?"p, - 2d\y,^cj, - [g^, + (e^^ + 1)A^A,] ^cj, + 2{e't' + l)A'^ A^^V 

+ (e^* - l)A"V(pA,) + (e^* - l)A(^F,f (31) 

The above transformations can be useful in writing the field equations in purely matter-frame form, or Einstein- 
frame form, although the result could be very complicated. One should not however use these transformations to 
write the actions for the three gravitational fields in matter-frame form as the transformations were derived under 
the assumption that the vector field is unit-timelike which is not valid at the level of the action. 

III. ROBERTSON- WALKER COSMOLOGY 

In this section I consider the evolution of the background cosmology. The field equations are first found in covariant 
form, to facilitate an easier transition to the inclusion of linear perturbations. A special coordinate system, common 
to calculations in cosmology, particularly in cosmological perturbations is introduced at the end of the section. 

A. Preliminaries 

1. Covariant description of Robertson- Walker geometry : geodesic congruences, metrics and projectors 

The background spacetime is assumed to be homogeneous and isotropic meaning that both of the metrics are of 
Robertson- Walker form. This assumption permits one to construct a smooth congruence of timelike geodesies, which 
are normal to a hypersurface of spatial homogeneity and isotropy. The words timelike and geodesic imply a metric and 
a compatible connection. Even though in this theory there are two different metrics (and compatible connections), 
this is not a problem. The homogeneity and isotropy is a property of the manifold and not the metric, and any of the 
two metrics may be used to construct such a congruence. 

Let the pair (7^,^, d^,) be the Robertson- Walker metric and associated metric-compatible connection in the Einstein- 
Hilbert frame. One can then identify the vector field A'^ as the unit vector field tangent to the congruence of 
timelike geodesies mentioned above, for the pair (7^j/,9p) (see appendix-IXI) . Since by construction A'^ is hypersurface 
orthogonal, the Frobenious theorem , geodesic equation and unit-timclike condition give 

= 0- (32) 

Now the Robertson- Walker metric is conformal to the metric of a static spacetime foliated by spaces of constant 
curvature of radius Vc which will be called the conformal static metric. Minkoswki space (which has an infinite radius 
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of curvature) and Einstein-static space (with positive radius of curvature) are two cases of conformal static metrics. 
This means that one can write "/^^ as 

Jfiu = b'^Vtii', (33) 

where b is the scale factor in the Einstein-Hilbert frame and where ry^j^ is the conformal static metric mentioned above 
with its associated metric-compatible connection. Similarly the matter frame Robertson- Walker metric with 
associated metric-compatible connection 9^, can be written as 

Itiu ^ a^Viiy^ (34) 

where a is the scale factor in the matter frame and 77^^ is also a conformal static metric with being its metric- 
compatible connection. 

Consider now a different unit-timelike geodesic congruence of curves with respect to {r]^^,dfj,). Let A'^ be the 
tangent vector field of this congruence assumed to be Killing (this can be accommodated by the symmetries of the 
static spacetime), which by construction obeys 

Vf.uAf'A'' = -1, (35) 

and 

A'^d.A^ = (36) 

where = rj^^A'^ . A further property of is that it is covariantly constant 

d^,A, = (37) 

which follows from the fact that it is Killing, geodesic and hypersurface orthogonal. Transforming 77^^ to 7^^ one 
finds that A^' is related to A^' by 

= ie-^A'', (38) 
a 

and 

A,,=ae-^A,,. (39) 

With the help of A^ or equivalently A^ one can construct two projection tensors, given by s^^ = —A^A^, which 
projects tensors along and g'^^ = 5'^^, -I- A'^Ai, which projects tensors on the hypersurface of homogeneity and 
isotropy. The two projectors have the property that s'^^A'^ = A^, q^a'fv — l^v^ ^-^id Q^v^^ — 0- Of cource in the 
case of a FLRW background q^^^ ~ q'^^ where 

q'^, = S'',, + A^A, (40) 

whereas when perturbations are included, q'^^ will aquire a perturbed part due to the perturbations coming from A^^ 
(see perturbation section) while q'^^^ is by definition unperturbed, and the two projectors will not be equal. 



2. Relations between the scale factors and between the conformal static metrics 

A further property of the FLRW spacetime is that any gradient of a scalar function will be in the direction A^ 
or equivalently A^^. Thus letting the background value of the scalar field be 0, its gradient is df^cj) — ~ ("a"^) ~ 

- (A''^^) A^. The same holds for any other scalar field, e.g. the scale factors in the two frames, or the background 

density and pressure of the fluids etc. 

The spatial metric in either frame is obtained using the spatial projector (|40|l acting on the corresponding metric. 
Thus in the Einstein-Hilbert frame the spatial metric is q^^j/ = 7^,^ + Af^Ai, whereas in the matter frame it is = 
+ e^'^A^A^. Equation then implies that the two are conformally related as q^i, = e~'^'^q^,y, which prompts the 
relation 



a = be-^. 



(41) 
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between the two scale factors. 

Using H33|l . (|34|l and (|39l) the relation between the two conformal static metrics and is obtained as 



with inverse 



(42) 
(43) 



3. Connections 

Let C^^j, be the connection tensor for the connection transformation 9^ — > 9^. Using (|29|l one gets 

C"^, = [-2<5"(^^) + e^^ry^^.A"] £ln& 

Let also the connection tensor which performs the transformation 9^ ^ 9^. Using (|29|l one gets 



COL 



-25", ..A 



£lna. 

A 



(44) 



(45) 



The reader is advised not to interprete as a partial derivative. It is a connection which annihilates the metric 7^^, 
i.e. da^fiv = without any reference to a coordinate system. 

Finally consider the connection tensor E°'^y for the connection transformation — s- d^. This is the same as 
dfi ^ dfi ^ ^ df^, and therefore E"^^ = —C"^^ + D"^^ + C"^^, where D"^^ is the connection tensor in 
adapted to the Robertson- Walker geometry. Alternatively E"^^ is obtained directly from (|29(l . I|42|l . (|43|l and l|37|l . 
The final expression is 



= -2 (^'^(9;,^ A'^A^A, 



(46) 



B. The field equations 



Here I consider the field equations adapted to the symmetries of the FLRW spacetime. I do not explicitely consider 
the vector field equation as it is trivially satisfied. 



1. The fluid tensors and evolution equations 

Lets first consider the fluid related variables, J, and S^y. The fluid velocity can be expressed in terms of A^ as 
Ufj, = aA^ and = which give A^u^^ = — e"^. For the same reasons as above (see apoendix- IXjl the fluid velocity 
is geodesic, i.e. u^^d^u^ = 0. The scalar source is then given by as 

J = e-^^{p + 3P) (47) 
where p and P are the FLRW background density and pressure of the fluid. The fluid current is given by (|23|1 as 

jf, = -2 sinh(20) e""^ a pA^, (48) 
and the fluid stress-energy tensor in the Einstcin-Hilbcrt frame as 

S^,u = a' [-(1 - 2e-^'^)pA^Ay + Pq^,] (49) 

where q^^^, = r]f_,xq\. 

Changing connection from 9^ to 9^, the energy conservation equation for the fluid (|26|l becomes 

A^9^p + 3p(l + w)A^9^1na = 0. (50) 

where w is the equation of state parameter such that P = wp. The momentum transfer equation H27|) is trivially 
satisfied. 
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2. The constraint equation 
Consider now the constraint (|20|) which when adapted to the symmetries of FLRW spacetime gives 

This equation is then inverted to get jl = /i(a, (p, A^^d^cj)). 



(51) 



3. The scalar field equation 

The vector field (|17|l adapted to the FLRW symmetries can be rewritten as — — f where f = A^F^ = 
2jl^{APdp(p). Using the above into the scalar field equation, gives a system of two first order equations, which are 



A'^df.r = -3{APdp In 6)F - 87rGae"^ J 



and 



APdp(f> = ^ae^'^F 

The two first order equations can then be combined into a single second-order equation which is 
1 



U 



3fl{A''d^,lnb){A''d^(l>) + AirGa^e-^^J + (^^^^.c^) (^A" d A^' dp In 



where 



Therefore to have a solution we need C/ ^ 0. 



4- The Lagrange multiplier 



(52) 



(53) 



(54) 



(55) 



Next in line comes the lagrange multiplier which is given by 

A = ^flV' - 167rG'sinh(2^)p 



(56) 



5. The generalized Einstein equations 
First lets compute the tensor Yp^. Using ifsB)) in ifTUIl one obtains 



Y,, 



3h(2</>)^U' - sinh(20)U + SttG (l - e"^^) p 



ApAf^ 



which when combined with H49|l gives the right hand side of the generalized Einstein equations as 



^e^^ ifiV -V) + SttGP 



+ 



sh{2(j>)flV' - sinh(2^)U + SvrG (e^'^'^p + p) 



ApAiy 



(57) 



(58) 



The conformal relation of the Robertson- Walker metric to the conformal static metric makes it convenient to use 
conformal transformations to calculate the Einstein tensor Gpi, for the Einstein-Hilbert frame metric (|33|l . Now the 
Ricci tensor of the conformal static metric is simply ^(jpu where K is an integer taking the value K = for a 
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spacially flat hypersurface (the conformal static space is Minkowski spacetime), K = 1 for a positively curved spatial 
hypersurface (the conformal static metric is Einstein-static spacetime) and K = —1 for a negatively curved spatial 
hypersurface. 

Then, after performing the conformal transformation the Einstein tensor is found to be 

^ he-^^A^A, - J - 2 [B^B, In b ~ {B^ In b) {B, In 6) - {r^Bc^Bp In b) t?^,] 



+ {B^\nb){Bf,\nb)fi^,. 
Changing connection to 9^ and using H42I) and (|43ll gives 
K 



Ze'^^A^A, - q^, +3 In 6 i^^+2e^^g^. 



B'^liib-iAf^Bplnb] [ 2A°'B, 



Giiv 



Contracting (|60|l and H58|) with A'^A'^ gives the generalized Friedmann equation as 



1 



A^dcAi-ib 



3(A^9^1n6)^ 



2 \ r^a^ 



while contracting the same equations with gives the generalized Raychandhuri equation as 



2d^ In 6 - 4:A°'Ba(l)A'^Bf3 In 6 - (A"5a In b)^ 
Other contractions give trivially zero. 



1 



e-^f'inV' -V)+e 



-40 



BttGP 



K 



(59) 



(60) 



(61) 



(62) 



C. Choosing a coordinate system 

A most convenient coordinate system that is commonly used in cosmological perturbation theory is the conformal 
synchronous coordinate system with t denoting conformal time and cc" the spatial coordinates. There are two choices 
regarding the frame for which the conformal static metric used above, takes the standard form. Since the connection 
to matter observables is through matter frame variables, it is more convenient to use the matter frame conformal 
static metric as taking the standard form, i.e. r/pp = — 1, ryp^ — and ?7.j, — q^^. This gives the matter frame metric 
as 



-dt^ + q^^idx^dx'' 



(63) 



in this coordinate system. The vanishing of the Lie derivative with respect to all the Killing vectors of the background 
spacetime gives 4> = 4>if) only, as well as Aa — Q and A°- = 0. The A^ component is found using H35|l which gives 
Aq = 1 and A^ = —1. Finally the components of fj^^ are fy^Q — — e""*"^, r/g- = and fj^i = 

Adopting the covariant equations of the previous subsection in this coordinate system gives the constraint equation 



as 



—a e 
2 



which must be inverted to get ^(a, ( 



Similarly the scalar field equation becomes 



1 

Tj 



b 



the Friedmann equation gives 



ie-^"? {^lV' + V) + e-^^ (sttGp ^ 



3K 



and the Raychandhuri equation 



Finally the fluid evolves according to 



^e'^'^ipV -V) + e-^'^ (87rGP+ ^ 



(64) 



(65) 



(66) 



(67) 



p + 3-(l + w)p = 0. 
a 



(68) 
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D. Effective Friedman equation and relative densities 

The physical Hubble parameter is as usual H = and after transforming the scalar field time derivative as 
(j) = Q rffna ^' effective Friedmann equation reads 

m' = 8nG,ff{p^ + p-^^^) (69) 

where the effective gravitational coupling strength is 

Geff = G -— (70) 

/-, I d4> \2 
^ dlna) 

and the scalar field density is 

Note that the effective gravitational strength is time varying. The relative densities fli for fluid i are then defined as 

1], = 87rGe//-^-— ^ (72) 
The above relation can also be used to define the relative density for the scalar field, O^. 



IV. PERTURBATION THEORY 

Cosmological perturbation theory dates back to the work of Lifshitz [ssf , who used a coordinate based approach and 
worked with the synchronous gauge. Many subsequent studies also adopted the same approach [s^. The synchronous 
gauge was found to contain spurious gauge modes |37| | causing confusion in some earlier studies as to what was 
the physical growing mode. Indeed some early studies identified these residual gauge modes and had to carefully 
remove them from the solutions. The existence of residual gauge freedom in the synchronous gauge lead Gerlach and 
Sengupta, Bardeen and others '38] to construct gauge invariant variables which as the name implies, were devoid of 
unphysical gauge modes. Covariant 39] studies of perturbation theory were initiated by Hawking and later developed 
by Ellis and Bruni and others into a fully covariant and gauge invariant theory. Gauge-ready approaches, were one 
can always choose a gauge at will, depending on what is more appropriate numerically were also studied by Hwang 
and Noh |4Qj. 

A. Perturbations of the gravitational variables 

1. Scalar field perturbation 

The scalar field is perturbed as 

(f>^4> + V (73) 

where ip is the scalar field perturbation. 



2. Metric perturbations 



The Einstein-Hilbert frame metric is perturbed as 

g^iK [ri^iu + hi_,^ (74) 
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where h^i, is the Einstein-Hilbert frame metric perturbation. The inverse metric is given by gf^" = ^ (f/^'^ — h^"^ 
where h'^" = fj^^ff^hap. One changes connection from to via V ^u,j = B^Ui, — (c^^ + f^J^ u\ for some form 
u^, where the connection tensor f^^, is given by 

/mi^ = ^V^^ {diJ-Vup + duffi^p - dpf]p,u) (75) 

Similarly the matter frame metric is perturbed as 

9tiv = (t?^^ + h^^) (76) 

where h^v is the matter frame metric perturbation. The inverse metric is given by g^^ = {t]^'^ — h^") where 
j^tiv _ ^iia^vp^^^_ Connection change from to on any form is as V^u^ = d^u„ — {c^^, + f^J^ u\ where the 
connection tensor /^^ is given by 

= Iv^" {d^V.p + d.v^p - dpv^..) (77) 

3. Vector field perturbations 
Let the vector field perturbation be a^, defined by perturbing the vector field as 

Af, = ae-* {A^ + a^) (78) 

and 

A^' = -e^ (A^^ + 0^") (79) 

where 

a" = (rj^'^a^ - h^""!^^ . (80) 
The field strength tensor F^^ then takes the form 

F^, = 2ae-^ 5[^a,] + A" (^0 - 4 In a) i[^a,] (81) 
Now define the "electric" and "magnetic" field parts of F^^ as and B^v which are given by 

= ^A^F^p = ^q'-^A^F^ff (82) 

and 

B^. = q"^q^.F„0 = -B,^. (83) 

The two tensors E^ and B^i, obey A^E^ — A^B^^, — 0, meaning that they are purely spatial. 
Explicitely they are given by 



d^iAf^ap) ~ Af^dpa^, + 1" (dc,4> - 9^ In a) (a^ + A^A^ap) (84) 
and 

B^, = 2ae-^r^fAc<^0] (85) 
which means that Ffjt^ = ae^'^ {A^Ei, — E^A^) + B^^. 
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4- Perturbation of the timelike vector constraint 

Lets condider the time-like constraint on A^^ which must be preserved even after the metric and the vector field 
are perturbed. This gives various relations between the metric perturbations and the vector field perturbations. 
Perturbing the constraint gives 

g^.A^'A'' = e*^ (r)^. + h^,) {A'' + {A" + a") (86) 

Expanding the right hand side, transforming Ty^j^ to ry^j^ with (|42|l and then imposing the constraint with l|19|l and 
gives 

h^.A'^A'' = -26-^"^ A^a^ (87) 

Two more relations are 

Af'af, = -Af.a^' (88) 

and 

h^-'A^A, = e^'^.A^A" (89) 

5. Relating Einstein- Hilhert and matter frame metric perturbations 
Now lets find a relation between h^^, and h^^. Start by perturbing the metric transformation which gives 

= 6-2"^%^ - 2 \e-^'^g^^ + 2 cosh(20) A^^^] - 4 sinh(2^)^(^M^) (90) 
Using H74|l . (|76|l . (|42|l and H78|l gives the required relation which is 

h^.^, = hf,^ + 2(1 - e-^^)A(^a^) + 2 [g^^ + A^^A^] ip. (91) 

while the inverse relation is 

= h'"' + 2{e*^ - l)A^''a''^ + 2 [q^"" + A'' A"] ip. (92) 

B. Perturbations of the fluid 

1. Perturbations of the fluid tensors 

The fluid density and pressure are perturbed a,s p — p + Sp and P — P + dP respectively. The fluid velocity is 
perturbed as 

u^^a {A^ + 0^) (93) 

and 

u'' = - (A^ + r ) (94) 
a 

with 

61^ = ry^-^ei^ - h'^'A^ (95) 
The shear is already a perturbation and obeys the identities Ty^'^S^^ = and ^^S^^ = 0. 
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2. Perturbing the timelike constraint on the fluid velocity 
This is similar to the timeUke constraint of the vector field, ft gives 

and 

A^'0^ = -A^Qf' 



(96) 



(97) 



Using (|91|l . one finds a relation between the fluid velocity perturbation and the scalar and vector field perturbations 



as 



which also gives 



Af^Of, = A^a^. - cp 
A'^u^ = -e^(l + if) 



(98) 
(99) 



3. Perturbing the fluid source tensors 

Using the above relations the scalar source perturbation SJ is found to be 

dJ = e~20 \^Sp + 3SP - 2{p + 3P)ip] 

Perturbing the current gives 

Sj^ = ae-^l - 2 sinh(20) [{p + P)d^ - Pa^ + 5pA^] 
-{e^^ + 3e-^^)p + 2 sinh(2<^)pj Af,ip 
Finally, perturbing the Einstein-frame stress-energy tensor yields 

6Sf,, = a^jpV -2Pg,,,(^- 2(1 -e-'*'^)pA(^a,) -2 [(l + 3e-4^)p + e-4'^F 
- (1 - 2e-^'^)5pA^A, + SPq^, + 2e-4^p + P)A(^0,) + E^, 

where S^^ = ^S^i.. 



Afj.A^(p 



(100) 



(101) 



(102) 



C. Perturbed field equations 

1. The perturbed fluid equations of motion 

Let the density contrast and sound speed be given as usual hy 6 = ^ and Cs — respectively. Then the perturbed 
fluid equations become 



and 



A^'^^6 + SiAf'd^ In a) (C^ - w)S + {1 + w) [ 9^0^ + ) = 



APdpe, + (1 - 3w)iAPdplna)e, - -d,{A'^A^Kp) 



(103) 



, Cl5 \ ^ APdpW, 



1 



l + w J l + w p + P 



= 



(104) 
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2. Perturbed Constraint equation 



The constraint equation yields 



(105) 



Unlike the unperturbed case, this need not be inverted. Rather it gives bp, directly in terms of the other variables 
which is then used in the relevant places in the other perturbed equations. 



3. Perturbed scalar field equation 



Following the same approach as for the background, start from (|17|l . which is then split into the background 
part f ^ found in a previous section, and a perturbation 7^, as = f ^ + 7'^. Then one performs the projection 
F'' = —TA^^ + q^^V^ where F = A^T^ . The scalar field F is again split into background and perturbed part as 

F = f + 7 where 7 = ae^^{a^f^' + ^^7'') is the perturbed part which makes 7** = (7 + TA^ay) + q^'^Y ■ 
After some calculations one gets the perturbation 7^ to be 



r 



^e'^iA'dx^) 



A-'h^',, - A^a^^A^" - 



o? V" 

while its projection on the hypersurface is 



{APdp$)A^ax + Af'dp^ 



A" 



L . 



The scalar field equation is then split into two first order equations given by 



-'S{APdpliib)j 



'^^Aip + F 



^A^dx V) - 87rGae-3^p [(1 + 3C^^)S - (1 + 3w) {A" a, + 2^)] 



where A = q'^'^dndi, is the Laplace-Beltrami operator and 



A''d,.ip 



1 

2U 



(106) 



(107) 



(108) 



(109) 



4- Perturbed vector field equation 



The divergence of F'^^ in terms of Ef^ and B^i, is 



V,,i^^, = —A^'d^E, 

a a 



r'^A^d^E^ + ^<f^LB^, + — (A^9^ ln6)£;. 



62- ■ a 

Its easier to perturb the vector field equation (|15() which does not include the lagrange multiplier. It gives 



(110) 



Kb 



+ e'^^r^df^B,^ + (i^4 ln6)^„l + q\\ n(A^dp4>) \d,ip + {Af'dp4))a, 



+8^Ga"(l - e-^^){p + P) {9, ~ a,) }=Q 



(111) 



which is a first order equation for E^. The other equation needed is a rearrangement of the definition of i?^ as a first 
order equation for a„ 



Af'd^ir^a,) = -E^ + q\ [dM^ap) + A^ {d^4> - 8^ In a) a,} 



(112) 
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while contracting with q'^al'^B gives 



2 

ttGo^P 



e-- {flV' ~V) + SttGP 



1 al pl^tJ-v 



{Cl5 ~ 2wip)qai3 + qf^^q^ptf,^ 



20 



qap 



(120) 



The above equation can be further simphfied by separating it into trace and traceless parts. The trace part is found 
by contracting with q'^'^ and is given by 



,20 



-e-- {flV' ~V) + SttGP 



q'^'h^^ + 247rGa2p(Cf (5 - 2wip) 



+3^e 



20 



a^V'APap + 2e^^{APdp4>)iA^dxip) 



(121) 



while the traceless part is 



-a —V —LLV — 



[(5 V + 87rG(55^,] = 



—Ll —V —LLV — 



-e^'t' {jiV -V)+ SirGP 



h 



(122) 



Now lets turn to the left hand side of the generalized Einstein equations. The perturbed Ricci tensor of ry^i/ + h^i, 



is simply given by 2d[\f\ , which gives the perturbed Einstein tensor as 



K 



and after expanding the connection tensors 
1 



(123) 



(124) 



The perturbed Einstein tensor in the Einstcin-Hilbert frame, is then obtained via a conformal transformation as 



SG^u = SHpu + 2 f'j^^da In 6 - ry^. 



2h'"^dadi3 In b + 2fj°''^ f^ijdx In b + h""^ A^Ap{APdp In 6) 
29^ Inb-e'^'^iA'^dplnbf 

Changing connection to (see appendix-Q and combining terms gives 

-^^M- = ^ [^''^hapv^^u - ihpj\ - - \r^dJih^, + i [ry^^^paA/i^^ - r'f^dj^, 

+ r^dJ(X)P + 2e"^ [d^ - &{APdp4>f] A'^AfKpApA, - e^^{APdp4>)A^dx [<tph^o) '7m- 
+ {APdp4>)\W^A^do.hf3x-APdih''^ 



(125) 



ApAjj 



+ e 



40 



A^aA/i^. - 2i"5(^/i,)„ + 2A''APA(^^d,)Kp 



26^^ I (^"9,^ 



+ e^^iA^dx In 6) lA^dihp, ~ 2AI'd^ Ji,)p ~ A{APdp4>)A^ A^KpApA, 



^^{APdp In 6) ^^9A/i"a - 2APr^dJifjp + ^{APdp$)Ao,APh''p 



„40 



2^2 In 6 - 4(A^aA In6)(^''9p0) - [A^dx l^bf] \hp, + A^A^'h'' pf,p. 



(126) 
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Now lets perform the contractions like in all the above cases. Contracting H126() with A" gives 



.-40 



e-^fq^^h,,, + SA^A^hf,, 
+ {APdplnb) [A^dx - 2APq''l'd^hi3p 

Contracting p26|) with A'^q'^^ gives 

-^A'-r Jm-^ + ^A^^q^qP^dpd^Kp - ^A>^q-J^d, + (qP^^^h^p 

^A^^q\Ah^, - {A^dx\nb)q\d, [A^A^^h'^^j) 
2d^ In b - 4(A^dx In b) (APdp^) - {A^dx In bf] A^f^hi"^ 



(127) 



+ 



(128) 



and likewise contracting (|126|l with q^a'f a yields 



ro.<tp5G^, = q^q'p^ 



q''^hpxqf^i, ~ itipy 



3ft-,, 



q\j''p)q''^dpdpKx - e'^q^^J''p^APA%dpKx 



Iro^Q^fAdr^h^ ~ ig^^r^jAft^, + ^e'^q'^^q''pAPA%dxhf.. + \ 



e^^APA^dpdx - r'r^M-hpX + 2e^^A-'dJp (q^'A^x) 

e^'^{APdp4>) [2A^qP''dphx. - A^dx (g^" V)] 



q\<tpAPdphp, - 2q\j^p.^APdpKp 



e^'^{APdp\ub) A^dxh^^ - 2APq^''dpKp + 2APdp [ApA'^h^ 



252 In 6 - '^{A^dx In b){APdp4>) 



- (A^dxlTibf [qP^ctphp. + ApA-'hP^q^p 
Taking the trace of the above equation, by contracting with cp^ yields 

q'^'SGp, = -^qP-'qP^dpdJl^x + 2e^^A''d,dp (qP^A^llpx) + V) - A (ApA'^hP, 



(129) 



^'^APA^dpdx [r^hp.) + 2e^HAPdp4>) [2A^r''dphx,. - A^dx [r^h 

^'^{APdplnb) [2A^dx {^qP-'hp^^+SAPdp (ApA'^hP^^-AAPqP'^d, 

hd^ Info - A{A^dx \nb){APdp$) - [A^dx \nbf] W^hpu + iApA^hP^ 



(130) 



which gives the traceless part as 
1 „„_ 



9 a9 /3- ga, 



5Gp, = 



- — ft^. + q'^dpdpKx - e^^APA^dpdpKx 



2^hp^ 



^dpd, [qP%x) + Idpd, {ApA^hP^) + ^e^'i'APA^dpdxhp 



„4</. 



„4<?i 



APdp<j> + APdplnb 



A^dxhp, - 2A^dph,x 



2d^\nb - A{A^dxliib){APdp4>) - {A^dxlnbf 



(131) 



The above contractions are then combined with their counterparts coming from the right hand side of the generalized 
Einstein equations. 
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V. IRREDUCIBLE DECOMPOSITION OF METRIC PERTURBATIONS 
A. Harmonic mode decomposition 

1. Einstein metric decomposition 

Lets write the Einstein metric perturbation into irreducible parts. This yields 

V = 2e-4'?SA^^,-2A(^r.)aaC-2A(^/i,) 

+ Ixq^. + (^r^g^. - Iq^.r^^ djpv + 2dJ(^q\) + Xi^u (132) 

The variables above are classified as follows: Scalar modes (5,C,x and v), vector modes [h^ and f^) obeying 
q^'^djiy = and g^^S^/i. = and tensor modes X/^i/ obeying q^^'Xi^.f = and q^'^dxXfi.iy = 0. 

2. Vector field decomposition 

The vector field perturbation is decomposed as 

= -EAf, + q\d,a + (3^ (133) 

with A'^fS^ — 0. It contains a scalar mode a, given by Aa — q'^'''^df^a„ and two vector modes /3p obeying ^'"'9^/3^ = 0. 
The "electric field" is also decomposed as 

= r^.^,E + (134) 
with E being a scalar mode given by AE — q^^d^E^ and two vector modes obeying as usual qf^^d^ei, — 0. 

3. Matter metric decomposition 

The matter frame metric is decomposed in a similar way as the Einstein frame metri as 

V - 2SA^A, - 2A(^g",)a„C - 2%/i,) 

+ lxq^.. + (r^g^. - ^9^-9"'') dadp,^ + 2aa/(^r.) + X^. (135) 

where A'^/i^ = = i'^x^'^ = 0. 

The variables above are classified in the same way as in the Einstein frame as follows: Scalar modes (5, C,x and 
v), vector modes (/i^ and J^) obeying qf^^d^h^ = and q^^d^f^ = and tensor modes Xt^v obeying q^^Xfj.f — and 
q^^dxx^. = 0. 

4- Relations bewteen different frame variables 

The relations between Einstein and matter frame modes can be read of from p32|) and (|135|l with the help of H91|l . 
They are as follows : 

(136) 

-4"^) a, (137) 

(138) 
(139) 

(140) 
(141) 
(142) 









c 


= c- 


(- 


X 


= x + 


6(p, 


V 


— V 






= - 


-(1 




— /aii 






— XfJ-v 
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5. Fluid velocity field decomposition 

The vector field perturbation is decomposed as 

= -EA^ + q\d,e + (143) 

with A^'v^, = 0. 

It contains a scalar mode : 6, given by A9 — q^^b^^Oy and two vector modes obeying q'^^^df^Vi, = 0. 

6. Fluid shear decomposition 



The shear perturbation is written as 

f:^,^{p + P)q%q^, 



(144) 



with A^afj, = A^a^i, = 0. The variable E is a scalar mode, the vector field a vector mode obeying g^'^o'p.iy = 0, 
while the tensor a^y is a tensor mode obeying as usual q^^a^^ = and q'^''(Jafj.,i, = 0. 

B. Gauge non-fixed equations for tiie scalar modes 

All scalar modes can be decomposed in terms of a complete set of eigenmodes of the Laplace-Beltrami operator. 
For example a variable A can be written as A{x'^) — J <Pk Y{x°',k-^) A{k-^), where the eigenmodes Y{x°',kf^) obey 
(A + fc^) y = 0. In the special case of a flat hypersurface with trivial topology, the eigenmodes are simply given 
Y — e*'^"^ and the integral transform above is a Fourier transform. The wavenumber k takes values depending on 
the geometry and topology of the spatial hypersurface. In the case of trivial topology, k takes values k — \Jk1 — 

where fc* is continuous obeying fc.^ > for a flat or negatively curved spatial hypersurface, and k^, — ^ where N is 
an integer obeying > 3 for a positively curved spatial hypersurface. Let us also choose the same coorcfinate system 
defined in section- ITTll 

1. Fluid equations 
The density constrast equation for scalar modes is 

5 = -3-(C2 - w)S +{l + w) i-k^e -\x + ) (145) 



a 

while the momentum divergence equation is 



.s - ^(1 - iyj)e + -^6 - -^e - ^ ( _ 3^ 1 2. (146) 



2. Scalar field equation 



The two equations equivalent to the scalar field equation are 

7 = -3^-f+^e'^^k^(ip + $a)+—fl$ 
a V J a 



-30; 



{1 + 3C^)S - {1 + 3w) (E + 2ip 



(147) 



and 



05 = ae ^7 - 



(148) 
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3. Vector field equation 

The scalar mode of the perturbed vector field evolve according to the two first order equations 
KbIe+^e\= -fi${ip - 4>a) + 87rGa2(l - e-^^)(p + P){e - a) 



and 



a = E-'E + 



a 



(149) 



(150) 



4- Generalized Einstein equations 
The scalar modes of the perturbed Generalized Einstein equations yield the Hamiltonian constraint equation 



lie 



(X + k'^v) + e 
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ae^^^ - Ksk'^E = SirGa'^p [5 - 2ip] 



the momentum constraint equation 

-^(^ + fc'^>) + ^ (2C + '>) - = 87rC?a2e-^^(p + P)^ + 2/i|<^ 
and the two propagation equations 

E = 2A-KGa^e-'^^p{Cl5 - 2wip) 



+ 2k' (e + e-'^ (k' - + ^'^) - 2 
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ic - 2fc2c] + 3^ae-^07 + 6 



J b' .jb 



and 



9 + e 



-i4> 



~ 1 1 9 

25- -y- -k 
3^ 3 







+ 2C + 2 





V + 2C 



= 167rGa2e-^^(p + P)S 



(151) 



(152) 



(153) 



(154) 



C. Gauge non-fixed equations for the vector modes 

Let m^, be an orthonormal triad of dual vector fields, normalized with respect to 77^,^ which give 

+ m^m^ + n^riy (155) 

and together with they form an orthonormal tetrad for the metric r]fxv 

Without loss of generality let be the direction of propagation of plane waves. Thus all vector modes are orthogonal 
to for example, — 0. Each vector mode X can then be decomposed into its two polarizations : 



X n,, 



(156) 



As it turns out there is no mixing between the two polarizations. Moreover they obey identical equations. The " +" 
and "— " labels can therefore be dropped without any confusion. 

Vector modes can also be decomposed in terms of a complete set of eigenmodes of the Laplace-Beltrami operator, 
just like scalar modes. The spectrum of the wavcnumber k is modified though, to reflect the spin one nature of the 

vector modes. In this case (again for trivial topology), k takes the values k = ^Jkl — 2^, where fc* is continuous 

obeying fc* > for a flat or negatively curved spatial hypersurface, and fc* = ^ where A'' is an integer obeying N > 3 
for a positively curved spatial hypersurface. 

Let us now find the equations for vector modes. 
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1. Fluid equation 



The vector mode fluid equation becomes 



a 1 + w 



2K 



(157) 



2. Vector field equation 
The two first order equations for the vector field are 



and 



(3 = e+[4>--]/3 



. b /,2 2K\ 
e+ -e+ A; + ^ e" 



6^/3 + 87rGa2(l - e-^^)(p + P)(u - /3) 



(158) 



(159) 



3. Generalized Einstein equations 



The vector mode momentum constraint is 



and the propagation equation is 



(160) 



(161) 



D. Equations for the tensor modes 

Using the orthonormal basis defined above, one can do a similar decomposition for the tensor modes into two polar- 
izations. The tensor mode perturbation decomposes into a basis which is written as symmetrized combinations 
of and n^. There are three possibilities, namely m^m,y, n^rii, and m^n^ + TOi^n^. However the traceless condition 
on the tensor modes, implies that the coefiicient of the first two must have opposite sign, hence there are only two 
independent polarizations given by 



(162) 



Since there is no mixing of polarizations again, the labeling can be dropped. 

Tensor modes are again decomposed in terms of a complete set of eigenmodes of the Laplace-Beltrami operator. 
The spectrum of the wavenumber k is again modified, to reflect the spin two nature of the tensor modes. In this 

case (again for trivial topology), k takes the values k = ^k1 — 3^, where fc* is continuous obeying A;* > for a flat 

or negatively curved spatial hypersurface, and fc* = ^ where N is an integer obeying A'' > 3 for a positively curved 
spatial hypersurface. 

The tensor modes then obey the equation 



H 




2K 



(T) 



(163) 
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VI. SUMMARY 



I have taken a covariant approach to formulate the Hnear perturbation theory about a spatiaUy homogeneous and 
isotropic spacetime. The covariant approach is particularly useful in theories with two metrics where there are two 
different metric compatible connections, one for each metric. 

The field equations were perturbed covariantly without adhering to a particular gauge or perturbation mode. This 
allows one to check explicitely that the equations are indeed invariant under infinitesimal gauge transformations. 
Mode decomposition was performed covariantly, and the equations for each perturbation mode were found, again 
without assuming a particular gauge. Special gauges for scalar modes are given in appendix While I have not 
considered the perturbed Boltzmann equation for thermalized fluids, this will remain unchanged when expressed in 
matter- frame variables. 

This completes the linear perturbation theory for Bekenstein's TeVeS theory about a FLRW cosmological back- 
ground. The equations presented here can be used to study the formation of linear structure and the Cosmic Microwave 
Background in this theory as was initiated in |2Q|| . 
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APPENDIX A: ANY UNIT TIMELIKE VECTOR FIELD IS GEODESIC IN AN FLRW UNIVERSE 

In this section I prove that that any unit timelike vector field in an FLRW universe obeys the geodesic equation. 

Let Qf^i, be the Robertson- Walker metric with scale factor a. Consider now a unit-timelike vector field tangent to 
a geodesic congruence of curves. As a property of Robertson- Walker metrics, is always orthogonal to a hypersurface 
of homogeneity and isotropy (see for example pi]). (Of course at least one such vector field exists, e.g. = V^t 
for some scalar function t e C°°M). Let x^, and be three unit-spacelike vector fields, which along with t^^ 
complete an orthonormal basis on TM. As a property of FLRW, they can be related to three linearly independent 
KiUing vectors of M as = ^^^ly = a^i2) ^^'^ ~ a^fs)' '^^^ above mean that 

t^'V^t^ = t^V^t^ = x^V^x^ = (Al) 

xA'V^x^ = (f ^'^") 

where in the last relation I have used the fact that Cf = for any function / <E C°°M. 

X 

Now the vector field 

A" = [l + c^f/'^tt' + cxf" (A3) 

is also unit-timelike by construction, for any choice of c e C°°M . The isotropy of Af , implies that there is no loss of 
generality in (|A3|) . 

Now consider A^V^A^. Using IjASp one gets 

^^V^A^ = c ci^ + (1 + c2)i/2a;^ (f 1^^) + + ^^)^''^ (^^V^x^ + x^V^t^) + {c\na^ U 
Now consider the term t^V^Xi/ -I- x^V^t^ which can be expanded as 

i^V^^a;^ + x'^V^.U = (fin a) + {y°'t^VpXc, + y^x^VfjQ y^ + {z"tN px^^ + z°'x^VpQ z^ 
where (|A1|) and (|A2|) have been used. However the coefficient of y^ in the above relation is zero, since 

y^t'^Vfjx^+y'^x^Vpt^ = y'^t^V px^ + y'' Ct^ ~ y^t^V ^xp 

= 2y"i^V[0X„] 
= 2y^/'dpe^ 



(A4) 
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The same holds for the coefficient of for the same reason and therefore 



ct^ + (1 + c^)^/^xJ £ln(ac) 



Therefore the choice c = ^ for any constant cq means that the unit-timehke vector field given by 



A^' = 



,2 1 1/2 



1 + ^ 



is geodesic. However any unit-timelike vector field can be related to by (|A6|) which completes the proof. 



(A5) 



(A6) 



APPENDIX B: GAUGE CHOICES 



1. Gauge transformations 



Consider a vector field generating a local one-parameter family of local diffeomorphisms (gauge transformations) 
Then under a gauge transformation, any tensor T transforms as 



T ^ T -I- rT 

where £T is the Lie derivative of T along 

Lets define a new vector field ^'^ by ^'^ = and — g^uSJ^ — a^/^- Now perform a split as 



(Bl) 



(B2) 



where A^^u 



0. The above vector field thus consists of two scalar modes ^ and '0 given by ^ = ^^Cm ^^i^ 



A?/; = q^'^d^^^i, and two vector modes which obey q'^'^d^uji, — 0. 

Now one can find the gauge transformations for all the perturbed variables. The scalar field perturbation transforms 

as 



1 



while the auxiliary scalar field perturbation 7 transforms as 



(B3) 



7' = 7 + 2^C/^ d''4>-{APdp4>? + {APdp4>){A''dAT^a) 



7 + 



6^ {A^'^^, In h) {A''d^(j3) + SttGc-^-? (p + 3P) 



The vector field perturbation transforms as 



which gives 



a' — a A — ^ 
a 

whereas (3^ is gauge invariant as expected. 

The vector field tensor i^^^ vanishes for the background meaning that _E, and B^^i, are all gauge invariant. 
The matter metric perturbation transforms as 



(B4) 
(B5) 

(B6) 
(B7) 



h' — h 



d(piu) + [APdplna) ^(^t) - {APdp In a) ^ 77^, 



(B8) 
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which gives 



X = 



h' 

r 



a 
2 

X- - 
a 

2 , 
V + -ip 
a 

1 
a 

/+- 
a 



A;V + 3(^''a^lna)^ 



h- 



A^'^^UJ-{A''^^\na)io 



where as the Einstein-frame metric perturbations transform as 

U> - 

^ 1 r 

C 



x' 



h' 
f 



- 1 
s + - 

a 
2 

X- - 
a 

- 2 , 
a 

+ i 

a 

UJ 

f+- 

a 



A'^d^.uj - {A''d^,\na)LU 



The Lie derivative of the fluid velocity is £u^ = and therefore the fluid velocity transforms as 
Now both the energy density and pressure are scalars given by p = u^u^T^" and P = ^q^^T^'^ and so 



5' = 

6' = 

--6P' = 

S' = 



S+ -{l+w){APdp\na)^ 
a 

p a 

S 

V 



3w(l + w){APdplna) - A^dpiv 



(B9) 
(BIO) 
(Bll) 
(B12) 

(B13) 
(B14) 

(B15) 
(B16) 
(B17) 
(B18) 
(B19) 
(B20) 

(B21) 

(B22) 
(B23) 

(B24) 

(B25) 
(B26) 



Using the gauge transformations above, a lengthly calculation shows that the gauge non-fixed equations derived in 
the previous section are all gauge invariant. This is a very powerful test that the equations are correct as given the 
complexity of the equations it is a very non- trivial matter. 



2. Conformal Newtonian gauge 
The Conformal Newtonian gauge is defined by 



X 

C 



-6$ 



(B27) 

(B28) 
(B29) 
(B30) 
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From the relation ^ = ^ + we also set S = and from x = X + S*^ we set x = —6$. Therefore the Einstein-Hilbert 
frame metric perturbations are given by 



c 



- ip 
^ - ip 

-(l-e-4^)a 



= 



(B31) 
(B32) 
(B33) 
(B34) 



a. Fluid equations 

The density constrast equation for scalar modes in the Conformal Newtonian gauge evolves as 

d = -{l + w) {k^e - 3$) - 2^{Cl - w)5 



where as the momentum divergence evolves as 



- I A;2 - 1 E + * 



(B35) 



(B36) 



h. Scalar field equation 

The two first order equations coming from the perturbed scalar field equation are 
.6 



7 = -3-7 + ^e-^^k"^ (^ip + ^aj +^fi^ -6^ - 2k'^C + SnGae-^^ Sp + 3SP +{p + 3P) - 2<pj (B37) 



and 



(B38) 



c. Vector field equation 
The scalar mode of the perturbed vector field equation is 

\ b 



Kb{e+-E\ =-fl(l){ip-(l)a)+87TGa'^{l-e-'^'''){p + P){0-a) 



and 



d = E + ^ + 



a 



(B39) 



(B40) 



d. Generalized Einstein equations 
The scalar modes of the perturbed generalized Einstein equations yield for the Hamiltonian constraint 



3$ + fc^C + 3^* 



+ - Ksk'^E = SttGo^p [5 - 2ip] 



the momentum constraint equation 



^+^C+r^-fi^^ = AnGa'^e-'^^ip + P)e 
ri 



(B41) 



(B42) 
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and the two propagation equations 



and 



+ +3^ae-^^-6 



b &2 



6l> + 3* + 2k'^C 



!> - § + e'''^ 



C+2 C 



(B43) 



(B44) 



3. Conformal synchronous gauge 

The conformal synchronous gauge is defined by S = and C = which fixes ^ = —{l — e~^'^)a and S : 
the standard notation lets also set x = h which gives x = h + 6(p and v = —^{h + 6r]). 



ip. Following 



a. Fluid equations 
The density constrast equation for scalar modes evolves as 

S = -3^(C2 - w)5 -{1+w) (k^0 + i/i 



while the momentum divergence evolves as 



.«(l_3«,)^+-^5--^. ^ 
' 1 + w 1 + w 3 



(B45) 



(B46) 



6. Scalar field equation 
The two equations equivalent to the scalar field equation are 

7 = -3^7 + ^e-^^k'^ U + ^o) + — 
a \ J a 

+ STrGae"^^ [5p + 2,5P - 3{p + 3P)ip] 



h + 6ip- 2k^Q 



and 



= ae -7 



(B47) 



(B48) 



c. Vector Bekenstein equation 
The scalar mode of the perturbed vector Bekenstein equation is given by the two first order equations 

Kb{e+^-^e\ =-fi${ip-$a)+87TGa'^{l-e-^^){p + P){0-a) 



and 



a = E- (p+[(j) ] a 



(B49) 



(B50) 
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d. Generalized Einstein equations 
The scalar modes of the perturbed generaUzed Einstein equations yield for the Hamiltonian constraint 



the momentum constraint equation 



a 



+ ae 



Ub 



J- Kek^E = 8TrGa^p[S -2<fi] (B51) 



2k^fi+^ 2k^C-h-6r] 

rf 



2A;2 ( ^ + /i^^ ^ + = ^TiGa'^e-^'^ip + P)k'^e 



and the two propagation equations 



-h-Q(p + 2k^C + ^e-^V + 2e-4^ ( k^ - ^) rj - 2j \h + 9ip - 2k^C 



-2(j) 



h + 6(p- 2kX 



,ib 



and 



h + Qi)- 2e-'^^k^'n - 2k'^C + 2 



b 

b^' 



h + 6f]- 2r c 



(B52) 



(B53) 



(B54) 



4. a-gauge 

In the a-gauge one sets the vector field perturbation to zero, a = 0, which fixes ^ = C- Since setting a = 
essentially fixes the gauge variable ^, it is no longer allowed to set S = 0. Therefore this gauge cannot be put in 
synchronous form. As a further gauge fixing condition, let u = 0, which can still be done, as the gauge variable ip was 
not fixed at that point. 



a. Fluid equations 
The density constrast equation for scalar modes evolves as 

6 = -3^(C2 - w)6 + {l + w) (-k^e -^x + k\ 



while the momentum divergence evolves as 



= -E-^{l-3w)9+-^6--^. , 



ri 



(B55) 



(B56) 



6. Scalar field equation 
The two equations equivalent to the scalar field equation are 

b a a 



+ SirGae-^^p [(1 + SC'^)S - (1 + Sw) {E + 2ip)] 



and 



(B57) 
(B58) 
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c. Vector field equation 

The scalar mode of the perturbed vector field equation is given by the first order equation 



i^s I i; + -£ = -fi^ip + 8TrGa\l - e-^^){p + P)e 



(B59) 



d. Generalized Einstein equations 
The scalar modes of the perturbed generalized Einstein equations yield for the Hamiltonian constraint 



3 V r. 



ae^'^cj)-/ - KbPE = SnGa^p [d - 2ip\ 



the momentum constraint equation 



-\x + - = 87rGa^e-^^{p + P)e + 2fi$cp 
6 ri 



and the two propagation equations 



-2 



E = 24wGa^e-'^^p{C'^6 - 2w(p) 



and 



-40 











+ 2C + 4 





C = 167rGa2e-^'^(p + P)E 



(B60) 



(B61) 



(B62) 



(B63) 



APPENDIX C: SPECIFICS OF THE PERTURBATIONS OF THE EINSTEIN TENSOR 

In this section I write out explicitely all the perturbed terms of the Einstein tensor to help the reader to follow the 
calculations. 

Changing connection to gives 

'A'^A^KfsA^A, - A^Aji^, 



+2(APdp4>) [r^A^do^hpxA^A, 



.40 



2A''A''dJip,A^ - A^dji^, + A'^A^d^hafiA,] } 



+Se^'l'{APdp^f [A^A^hx. - 2A''A^K,3A^A^ 

= d^d,h\ + 2{APdp^)A^A,A^dih° 



(Cl) 
(C2) 



d'^h^. = r^dJih^, - 2e^'^i^2d^A(^A^K)x + {A'dp^) A^'dxh^. - 4A"^^4ft/3(^^.) 
+4{APdp^f \A"A'^hc,0A^A, - 2A''A(Xw 



(C3) 
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d^4> - QAPdp4>f 



d''h\ = vP^dpdxh\ - 2iAPdp^)e''^A^dxh'-, 



(C4) 
(C5) 



2f^^dM = e^'i'iA^dxlnb) Ia^ Bph^, ~ 2Af' df^^^p ^ A{APdp4>)A'' A^KpA^A, 



(C6) 



2h°'^d^dp In 5 = 2e^^A°'A'^Kp [-9^ ln6 + 2[A^dx In 6) (^"9, 



(C7) 



2fi°'f^fKdxlnb = e'^'''iA''d^\nb) 



APdph\ - 2A^d^h"f, 



e-^^iAPOp In 6) A^dxh'^a - 2APf]°'P daJipp + ^e'^^{APdp4>)A°'A^Kp 



(C8) 
(C9) 



52 In 6 = e^"^ In 5 - 2{A^dx In 6) (A^^p 
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